On a principle of cosmological uncertainty 
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Abstract 

We show that cosmological observations are subject to an intrinsic uncertainty which can be 
expressed in the form of an uncertainty relation similar to the Heisenberg principle. This is a 
consequence of the fact that the four dimensional space-time metric information is projected into 
the one-dimensional observational red-shift space, implying a limit on the amount of information 
which can be extracted about the underlying geometry. Since multiple space-time configurations 
can lead to the same red-shift, there is an unavoidable uncertainty about the determination of the 
space-time geometry. This suggests the existence of a limit about of the amount of information that 
cosmological observations can reveal about our Universe that no experiment could ever overcame, 
conceptually similar to what happens in quantum mechanics. 
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Cosmology has some distinctive features respect to other branches of physics: we can 
observe the Universe only from one point, and we are unable to probe directly the geometry 
of space-time on large scales. We can only detect particles or radiation which have reached 
us, and extract from these observations information about the space-time between us and 
the sources. Most of cosmological scale distances can in fact only be deduced from red-shift 
measurements, assuming an underlying space-time metric. The Hubble law for example is 
based on the assumption of large scale isotropy and homogeneity of the metric describing 
the Universe. We know nevertheless that the Universe had a certain level of inhomogene- 
ity at very early times, as predicted by inflation, which is supposed to be the seed for the 



later process of structure formation due to gravity. There has been an long debate jl-5| 
on whether this perturbations of the homogeneous space-time metric could explain cosmo- 
logical observations, in particular if they could induce effects equivalent to a cosmological 
constant or dark energy. In this letter we consider what are the statistical effects of these 
inhomogeneities on the determination of background cosmological parameters, focusing on 
their complementary role rather than looking at them as alternatives to each other. Let's 
consider a Friedman universe with scalar perturbations. In longitudinal gauge the metric 
can be written as 

ds 2 = -(1 + 2^)d7] 2 + (1 - 2^)5 ij dx i dx j . (1) 

Dackground to be flat and the perturba- 



Sfl a photon is emitted from a source 



where we have assumed assumed for simplicity the 
tions to be of perfect fluid type. Lets us suppose [6 
located at a comoving coordinate r$ and time ts with a wave length As, and is reaching an 
observer at the center ro = 0, with a wavelength Ao • Its total redshift can be expressed 
as the sum of two components . One is due to the homogeneous background expansion, zh, 
and another due to the perturbations of the gravitation potential, zj : 

^ = i + z = i +Zjr + Z/ = ?M(i + [^v.n]°-2 ( X ° d\i,) , (2) 
As a{t S ) J\ s 

Z H = 77 r - 1 , 3 
a(t s ) 

Zl = (1 + z H )(ty + v ■ n]° s - 2 f X ° d\^ ) . (4) 

In the following we will neglect the contribution coming from the difference in the potential 
and the peculiar velocity between the observer and the source to get 

/■Ao * . 

Zl = {l + z H ){-2 d\4> ) w -2(l + z H )Ar$> , (5) 

Jx s i=l 



where in the last approximate equality we have chosen the comoving coordinate as the affine 
parameter along the null geodesic, we have divided the interval (rg, 0) into N subintervals 
of equal comoving length Ar and approximated the integral with a summation. 

The time evolution of the Fourier transform of the gravitational potential is given by: 

where D(rj) is the growth factor, T(k) is the transfer function, and (j)p(k) is the primordial 
perturbation. For a matter dominated universe the time derivative of the gravitational po- 
tential is zero, while for a Universe with cosmological constant there can be a significant 
contribution, which in the case of the cosmic microwave background radiation for exam- 
ple leads to the so called integrated ISW effect. While there is a redshift dependence for 
cf/ ( a^f) ' we w * u ne gl ec t it here, and assume an average constant value a along the geodesic. 
This implies that we can pull this constant factor out of the summation, and the only quan- 
tity summed is proportional to the primordial perturbation of the gravitational potential. 
Assuming this to be gaussian, we can approximate it as a discrete random variable ipi with 
uniform probability of having a positive, zero or negative value, corresponding to the stan- 
dard deviation of the gaussian field a^. This correspond to approximating the gaussian 
distribution with a uniform distribution with the three discrete values ipi = {<J^, 0, — a^}. 
Finally we to get 

N 

zj = -2{l + z H )aArJ2i>i , (?) 

i=l 

where ipi is a random value in the % — th sub-interval determined by the time independent 
part of the i()(k, rj), i.e. by ipp(k)T(k), and k is the momentum corresponding to the scale Ar. 
For large scales T(k) is approximately 1, so ipi is of the order of the primordial spectrum 
ipp{k). From the above assumptions we can derive the probability distribution of zj. In 
this discrete model every space-time configuration along the photon path corresponds to 
assigning a set of N values to the random variable fa. 

Let's call {n + ,n_}, the number of positive and negative values of the random variable 
fa. The corresponding Zj will be given by 

zi = -2(1 + z H )aAra 4 ,(n + - ra_) = Az min S , (8) 
Az min = -2(1 + z H )aAra^ , (9) 
S = (n + -n_), (10) 



where we have introduced the Az m i n , the minimum possible inhomogeneous contribution to 
the total redshift. According to our discrete model the total inhomogeneous redshift is a 
finite multiple of this minimum contribution, so we get: 



-N<^<N, (11) 

where the lower and upper bounds correspond to configurations of all positive or all negative 
values of ipi, i.e. n + = S or n_ = S. Since every interval can have three possible values of 
ipi, we have a total of 3^ possible different configurations. Because of the fact that what it 
matters is the sum (n + — n_), but not the order in which positive or negative perturbations 
appear, different space-time configurations can lead to the same value of Zi. This is the origin 
of the intrinsic uncertainty on determining the origin of the total redshift, and consequently 
to estimate other background cosmo logical parameters such the Q\. 

The number of space-time configurations corresponding to a given value of zi can be 
calculated analytically. Let's start from the number of configurations corresponding to a 
given set n + } of positive and negative values of ipi, which is is given by: 

Nc(n-,n+) = . — ■ v . (12) 

n^\n + l{i\ — — n + )\ 

The number of different configurations corresponding to a given value of zj, i.e. to a given 
value of S — (n+ — n_), can be computed in the case of even N and S > as 



(N+S)/2 

N s = ]T N c (n+-S,n+) = 2 F 1 

n+=S 



[l + S-N)/2; (S-N)/2;S+1;A 



s\(N-sy. 



(13) 



where 2 -Pi is the Gauss hypergeometric function. Since the distribution is symmetric in S 
respect to its central value, a similar expression can be found for S < 0. The formula derived 
above is consistent with the total number of possible configurations, i.e. we have 

E N S = 3 N . (14) 

S=-N 

The figure 1 shows a plot of the expected probability P(S) and P{zi) for different values 
of N. As it can be seen the dispersion around the central value (corresponding to zj — 0), 
is increasing as N increases, while its average value remain zero because of the symmetry 
of the distribution. Since this dispersion correspond to an uncertainty on the estimation of 
zh from the total redshift z because zh = z — zj, we deduce that an increased accuracy 
in estimating the effects of the inhomogeneities, i.e. a larger N, correspond to a larger 
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FIG. 1: The probability P(S) and P(zj/Az m , n ) are plotted for different values of N. The thick line 
is for N = 30, the dotted line is for N = 20, and the dot-dashed line is for N = 10. The difference 
between the two distribution is due to the fact that as N increases Ar decreases, lowering Az m i n . 
As it can be seen the dispersion increases with N. 

uncertainty on Zh, and consequently on the background cosmological parameters on which 
zh depend. Schematically we can write that 

Azj oc 1 , (15) 
Az H oc N , (16) 

(17) 

where the proportionality symbol has to be interpreted qualitatively, since we have not 
quantitatively estimated it, but we have inferred it heuristically from the dispersion of the 
expected probability distribution P(zi) and from the intuitive argument that a larger number 
of subintervals N corresponds to a more accurate modeling of the perturbations between the 
source and the observer, i.e. a smaller error Azj. We can express this with an uncertainty 
relation of the type: 

AzjAz H > f(r 8 - r ) , (18) 

where f(r) should be a monotonous increasing function of the comoving coordinate, since 
the uncertainty increases with the comoving distance of the source, due to a larger number 
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of possible space-time configurations. Alternatively we could for example write it as: 

A/AO A > f(r) , (19) 

where AI stands for uncertainty on the determination of the perturbative part of the metric, 
and Q\ is supposed to encode all the information about the homogeneous background, i.e. 
we are neglecting the effects of the curvature. A more accurate treatment could involve 
the use of stochastic calculus or of a random variable which can take more values following 
different physically motivated probability distributions. Nevertheless we expect that such 
improved approaches would not affect the the existence of such a lower bound, but could 
possibly modify its quantitative estimation. 

The dependence of the lower bound on the comoving distance of the source can be in- 
terpreted in terms of information theory, as the consequence of the fact that from a single 
piece of information, i.e. the ratio between the energy of the photons at the source and at 
the observer, we are trying to extract an increasingly larger amount of information, i.e. the 
space-time configuration between the source and the observer. Since the number of possible 
space-time configurations goes like 3^, if we keep Ar constant, we are dealing with an in- 
creasingly larger number of possible space-time configurations. Vice versa if we are keeping 
N constant and increase r$ we will be able to resolve space- time inhomogeneities only at 
larger comoving scale rg/N, so A J will increase. 
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